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Abstract 

We study the gross properties of even-even and N = Z nuclei with A = 
48 — 64 using shell-model Monte Carlo methods. Our calculations account for 
all Ohui configurations in the /p-shell and employ the modified Kuo-Brown 
interaction KB3. We find good agreement with data for masses and total 
B{E2) strengths, the latter employing effective charges = 1.35e and e„ = 
0.35e. The calculated total Gamow-Teller strengths agree consistently with 
the i3(GT+ )-values deduced from (n.p) data if the shell model results are 
renormalized by 0.64, as has already been established for sd-shell nuclei. The 
present calculations therefore suggest that this renormalization (i.e., Qa = ^ 
in the nuclear medium) is universal. 
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I. INTRODUCTION AND MOTIVATION 



The interacting shell model [1] is generally considered to be the most fundamental theory 
of the nucleus short of an explicit solution of the A-body prol )lom. Indeed, it is the conceptual 
basis for most other nuclear models. But in contrast to the atomic shell model, the residual 
interaction among the valence nucleons plays an important role. This interaction can be 
determined either by adjusting its matrix elements to fit a large volume of data, (as been 
done successfully for the p-shell [2] and sd shell [3]) or from the nucleon-nucleon G-matrix 
[4]. 

Due to the combinatorial increase of configurations with the numbers of valence nucleons 
and orbitals, the conventional solution of the shell model by matrix diagonalization has been 
limited in heavier nuclei to truncated calculations. For example, complete Ohoj calculations 
in the /p-shell have been possible only for nuclei with A < 48 [5] and, very recently, for the 
ground state energy and Gamow- Teller strength of the nucleus ^"^Cr [6]. 

Conventional studies of /p-shell nuclei with A > 50 must be performed in severely 
truncated model spaces. This is unfortunate as, for example, the nuclei in the iron region 
play a crucial role in a supernova, when the Gamow- Teller strengths determine the electron 
capture rate and thus the dynamics of the early collapse. Truncated shell model calculations, 
however, recover only in part the well established Gamow- Teller quenching [7-9] and are 
thus inadequate for predicting the Gamow- Teller strengths of those nuclei, for which no 
experimental information is yet available. 

Nevertheless, the conventional calculations performed for ^ < 50 support the conjecture 
that the shell model, if performed within a complete Ohcu basis, is able to describe the fp- 
shell nuclei. In particular, the systematic studies of the A = 48 nuclei [5] indicate that the 
residual KB3 interaction [10] is well suited to complete /p-shell calculations. This interaction 
has been derived by minimally modifying the monopole strength in the original Kuo-Brown 
interaction [4]. 

The results we present in this paper support that conjecture and demonstrate that com- 
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plete OfiLU calculations performed with the KB3 interaction are capable of describing the 
ground-state properties of even-even and N = Z nuclei in the /p-shell with A < 64 (for 
larger A the (7g/2-orbital must be included). Our studies make use of the recently developed 
shell model Monte Carlo (SMMC) method [11] in which few-body observables are calculated 
at finite temperature. The SMMC exploits the fact that most of the billions of configura- 
tions in these nuclei are quite unimportant for general nuclear properties, so that only a 
subset of the relevant configurations need be sampled. In comparison with the conventional 
diagonalization method, the SMMC is not yet capable of providing detailed spectroscopic 
information. However, it has been shown that the SMMC approach is well-suited for studies 
of both ground state properties [11-13] (obtained in the low-temperature limit) and nu- 
clei at finite temperature [14]. In this paper we systematically study important properties 
of /p-shell nuclei, including masses, and Gamow- Teller, E'2, and Ml total strengths. As 
our results demonstrate, these calculations are in agreement with experiment for all these 
quantities over the wide range of nuclei studied (A = 48 — 64) . 

II. THE SHELL MODEL MONTE CARLO METHOD 

The shell model Monte Carlo method is based on a statistical formulation. The canonical 
expectation value of an observable ^ at a given temperature T is given by [15,11,12,16] 
(/3 = 1/T) 

where U = exp{—(3H) is the imaginary-time many-body propagator and Tt^U is the canon- 
ical partition function for A nucleons. The shell model Hamiltonian H can be cast in the 
form 
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where are the single particle energies and Oa represent a set of one-body density operators 
(O denotes the time-reverse of O). The Hamiltonian in Eq. (2) is manifestly time-reversal 



invariant if the parameters Va that define the strength of the residual two-body interactions 
are re ah 

The key to the SMMC method is to rewrite the propagator [/ as a functional integral 
over one-body propagators. To achieve this goal, the exponent in U is split into Nt time 
slices of duration Af] = jS/Nt, 

U=[e-^^^f. (3) 

The many-body propagator at each time slice is linearized by a Hubbard-Stratonovich trans- 
formation [17]; i.e., it is transformed into an integral over a set of one-body propagators that 
correspond to non-interacting nucleons in fluctuating auxiliary fields defined by complex c- 
numbers cTq^ (n = 1, Nt). The expectation value of A then reads 

, . _ TT4Ae-f^^) ^ JD\a]W{a)<!>{a){A)^ 

^ TrA(e-/3^) ~ J D[a]W{a)<^>{a) ' ^ ^ 

where the metric is 

D[a] = n„„[da„„rf<„A/?|y„|/27r], (5) 
and the approximation becomes exact as A^j — )■ oo. The non-negative weight is 

W{a) = \C{a)\exp{-WK\\a^„\^A^), (6) 

an 

where ^(cr) = TtaU^ is the partition function of the one-body propagator = U^^ ■ ■ ■ Ui, 
with U„ = exp{—A(3h„) , and the one-body Hamiltonian for the n^^ time slice is 

K = Y^iC + SaVa<7c,n)Oa + (e« + V^CT^^) , (7) 

a. 

with = ±1 for Va <0 and s^, = ±t for > 0. The "sign" is = ^(cr)/|^(cr)| and the 

expectation value of A with respect to the auxiliary field a is 

{A). = TTAAUjC{a). (8) 

Both ^(cr) and {A)a can be evaluated in terms of the matrix that represents the evolution 
operator in the space of Ng single-particle states. In the applications discussed below the 



trace is canonical corresponding to a nucleus with a fixed number of nucleons [16]. Details 
of the transformation from the residual particle-particle interaction to the used above 
can be found in [11]. 

If all Va < 0, then the sign is ($) = 1. However, for realistic nuclear interactions such as 
KB3, about half of the V^'s are positive generating a sign problem (where the uncertainty in 
$ is larger than ($)). To overcome this problem, we extrapolate observables calculated for 
a family of good-sign Hamiltonians Hg (with g < 0) to the physical Hamiltonian at g = 1 
[12]. 

The results presented in this paper correspond to various observables of the nuclear 
ground state: the energy {H), the total B[M1) and B[E2) strengths, and the Gamow- 
Teller strength. The total B[M1) strength is defined as B[M1) = {fi~). where the magnetic 
moment Jl is given by /i = X]i I^n |<7/' + dsS^^i I^n is the nuclear magneton and gi^ g^ are the 
free gyromagnetic ratios for angular momentum and spin, respectively (gi = 1, gs = 5.586 
for protons, and gi = 0, Qs = —3.826 for neutrons). The total B{E2) strength is given 
by B{E2) = (Q^), where the quadrupole operator is defined hy Q = CpQp + e„Q„ with 
Qp{n) = ^p{n) Si ?'i^^2(^i, the sum ruus over all valence protons (neutrons). The effective 
charges ep^„ account for coupling to the giant quadrupole resonance outside our model space. 
For the oscillator length, we used b = I.OIA^^^ fm. The total Gamow- Teller strength is given 
by B{GT±) = (((TT-i-)^). We also explore the isovector monopole pairing in the ground states, 
as described below. 

The SMMC studies presented below have been performed in the complete set of 
0/7/2,5/2 — lp3/2,i/2 configurations using the modified Kuo-Brown KB3 residual interaction 
[10]. Each calculation involved 4000-5000 Monte Carlo samples at each of six values of the 
coupling constant g equally spaced between —1 and 0: extrapolation to the physical case 
(g = 1) was done by the method described in Ref. [14]. In many cases the assumption 
of a linear dependence on g resulted in acceptable fits (x^ per degree of freedom less than 
one). However, we have estimated the errors in our calculation conservatively by adopting a 
quadratic extrapolation in g. For the Hamiltonian we have made use of the variational prin- 
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ciple which ensures that (H) has a minimum at g = 1 [12]. The calculations were performed 
at (3 = 2 MeV~"^ (which is sufficient to guarantee cooling within a few hundred keV of the 
ground state for even-even nuclei [14]), and Nt = 64 time slices (which results in negligible 
discretization errors). 

III. RESULTS 
A. Comparison with direct diagonalization 

A few of the nuclei we consider here with SMMC methods have previously been studied by 
direct diagonalization using the same residual interaction. In all cases, the results obtained 
by these two very different methods agree. 

Poves and collaborators calculated Gamow- Teller strengths of B{GT^) = 1.263 for ^'^Ti, 
4.13 for 4^Cr [5], and 3.57 for ^^Cr [6], while our SMMC calculations yield 1.13 ± 0.18, 
4.37 ± 0.3-5 and 3.51 ± 0.27 for these nuclei, respectively. For the two iron isotopes ' 'Fe 
and ^^Fe, Caurier et al. [6] calculated the Gamow- Teller strengths in a series of direct 
diagonalizations with a decreasing level of truncation. In Fig. 1 we compare the results of 
these truncation series with our full Ofico SMMC results. Our results (6.05 ±0.45 for ^^Fe and 
3.99 ± 0.27 for ^^Fe) agree with the values obtained by Caurier et al. upon extrapolation 
to no truncation (5.5 ± 0.5 and 3.7 ± 0.5 for ^"'Fe and ^^Fe, respectively). Fig. 1 also 
clearly demonstrates that complete (Ohu) shell model calculations recover significantly more 
quenching of the Gamow- Teller strength than truncated calculations. Note that the single 
particle estimate for the total Gamow- Teller strength in both nuclei is B{GT^) = 10.3. 

For the total B(E'2) strength from the ground state we can compare our results with 
those of conventional shell model calculations for the nuclei "'^Ti [5] and "'^Cr [18]. Using the 
effective charges = 1.5e and e„ = 0.5e for protons and neutrons, respectively [5], these 
authors find B{E2) = 583.5 e^fm^ for ^^Ti, while we calculate 630 ± 55 e^fm^. For later 
reference it is noteworthy that roughly 80% of the strength lies in the transition from the 
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ground state to the first state [5], in agreement with experiment. For ^^Cr, Poves et al. 
calculate the total B{E2) strength of 1280 e^fm'' [18], in agreement with our SMMC result 
of 5(^2) = 1375 ± 110 eW. 

/.From Ref. [5] we deduce the energies (H) = -24.5 MeV and -32.9 MeV for ''^Ti 
and ^^Cr, respectively. Our SMMC results for these quantities are —23.9 ± 0.4 MeV and 
-32.3 ± 0.4 MeV, respectively. For ^°Cr we calculate {H) = -40.0 ± 0.4 MeV, while Caurier 
et al. [6] give —40.57 MeV. The slight underbinding we find is expected, as our calculations 
have been performed at the finite temperature T = 0.5 MeV and therefore should contain a 
small excitation energy of a few hundred keV [14]. 

More details of the comparison of SMMC calculations with exact diagonalization results 
with realistic interactions are given in Ref. [19]. 

B. Energies 

In this subsection we compare our calculated binding energies for the various nuclei with 
experimental data. The Coulomb energies, which are not included in the KB3 interaction, 
are calculated (in MeV) as [5] (tt = number of valence protons, p = number of valence 
neutrons, n = tt + p) 

i?Coui = ~ ■ 0.35 - 7TP ■ 0.05 + TT • 7.289. (9) 

Following Ref. [5] we have added to the calculated energy expectation values an energy shift 
of 0.028-n(?7 — 1) MeV to correct for a "tiny" residual monopole defect in the KB3 interaction. 
Note that the KB3 interaction uses the single particle energies from the experimental levels 
in ^-^Ca. Thus, our energy scale is set by e = —8.363 MeV for the fjf2 orbital, obtained from 
the neutron separation energy of ^^Ca. 

The nuclear mass excesses AM obtained this way (relative to "^^Ca) are compared to the 
data in Fig. 2a. The results obtained are very satisfactory, as can also be seen in Fig. 2b, 
which shows the discrepancy SAM between our calculation and data. Considering that our 
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finite-temperature calculation includes a small internal excitation energy of a few hundred 
keV, the reproduction of the mass data by the shell model using the KB3 interaction is 
better than is suggested by Fig. 2b. 

Two remarks can be made about our energy results. First, the accuracy of our choice 
of Coulomb energy (9) can be estimated by comparing the mass differences AM for the 
mirror nuclei (^^Fe, ^*^Cr), (^"^Fe, ^''Cr) and (^^Fe, ^^Ni). We find deviations between 0.1 
MeV (A = 48) and 0.5 MeV (A = 54) showing that our parametrization is sufficiently 
accurate, but might become less reliable with increasing Z. Second, in contrast to the other 
nuclei studied, the SMMC overbinds the Nickel isotopes, as well as ^^Fe and ^^Cr. As these 
nuclei have semi-magic proton or neutron numbers [N = 28), our results indicate that the 
KB3 interaction slightly overemphasizes the = 28 shell closures. Indeed, the overbinding 
is strongest for the double-magic nucleus "'^Ni. 

C. The B{M1) and B(E2) strengths 

We have calculated the total B{M1) strengths of the various nuclei using free nucleon 
g-factors, as listed in Table 1. Unfortunately the known Ml transitions in the /p-shell nuclei 
are mainly to low-lying states that exhaust only a small fraction of the total strength, so that 
a comparison to data is not possible in most cases. However, for a few of the nuclei studied 
here, Richter and collaborators [20-22] have used high-resolution electron scattering to study 
the B[M1) strength by means of high-resolution electron scattering in an energy window 
large enough to contain most of the strength. For the = 28 isotones (in units of /i^) they 
find B{M1)= 4.5 ± 0.5 (^°Ti), 8.1 ± 0.8 (^^Cr) and 6.6 ± 0.4 (^^Fe) for excitation energies 
between 7 and 12 MeV [21]. If one considers that this energy window should contain about 
75% of the total strength [21], our SMMC results (5(M1) = 12.5 ± 1.0 (^^Ti), 18.9 ± 2.2 
(^^Cr), 16.5 ± 2.8 (^"^Fe)) are roughly twice the observed B{M1) strength for these nuclei, 
supporting the idea that the spin g-factors are renormalized in the nuclear medium [23]. 

We note that both the independent particle and Ip — Ih shell models, discussed in Ref. 
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[21] predict a linear rise of the total B[M1) strength with the number of valence protons in 
the = 28 isotones. Our SMMC calculation does not support this trend, as the calculated 
B{M1) strength in ^''Fe is no larger than that for ^^Cr, in agreement with experiment 
[21]. We also note that the SMMC results might be tested by comparing the isotone pairs 
(^^Cr,^'*Fe) and (^^Cr,^^Fe) - the calculations predict roughly the same B{M1) strength in 
the Iron isotopes, but a significantly smaller B(M1) strength in ^"^Cr than in ^^Cr. 

For s^Ni Mettner et al. determined a total B{Ml) strength of 16.9i|i while the SMMC 
result is 20 ± 2. Within the large (experimental and theoretical) errors, these results are 
compatible with a significant renormalization of the (spin) g-factors. 

Table 1 also lists our calculated total B[E2) strength for transitions from the ground 
state. We have calculated these quantities by adopting the same effective charges used in 
a recent truncated shell model calculation of ' 'Fc [7] (e^ = 1.35e and e„ = 0.35e). For 
comparison Table 1 also lists the measured B(E2) values for the OJ*" — )■ 2i transitions in 
these nuclei; in even-even nuclei this transition typically exhausts about 70 — 80% of the 
total strength. For example, from the (e, e') data [28] we calculate ratios between the total 
observed B{E2) strength and the experimental B{E2,0'l — >• 2^) by factors of 1.35, 1.14 
and 1.20 for the Nickel isotopes ^*^'60.64^j^ respectively, while our SMMC results exceed the 
experimental B(E2, OJ*" — >• 2'f ) values by 1.38, 1.14 and 1.53. The overall level of agreement 
is illustrated in Fig. 3, where we compare the total calculated B{E2) strength with the 
experimenal B{E2, 0'^ 2'^) values for those /p-shell nuclei with either semi-magic proton 
or neutron number N = 28. A similar comparison has been presented in Ref. [29] based on 
a strongly truncated shell model calculation and a different residual interaction. To achieve 
overall agreement with the data, Ref. [29] used somewhat larger effective charges (e^ = 1.4e 
and e„ = 0.9e) in order to compensate for correlations missing in the truncated model space. 
Note that the use of the larger effective charges would increase our total B[E2) strengths 
for the nickel isotopes by more than 60%. 

While our calculation apparently compares nicely with data for the nuclei in the middle 
of the shell, our B{E2) values as listed in Table 1 appear too low for some Chromium 
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|-48,50Qj,^ and Titanium isotopes (^^Ti). For these nuclei larger effective charges are required, 
perhaps indicating a greater importance of sd-shell configurations at the beginning of the 
shell. 

D. Gamow-Teller strength 

In previous publications [12,13] we have shown that full Ohuj shell model calculations 
recover significantly more quenching in /p-shell nuclei than truncated 2p — 2h calculations; 
a finding that is in agreement with the recent work of the Strasburg-Madrid group [5,6]. 
Our work suggested that isoscalar proton-neutron correlations [30] and proton and neutron 
pairs coupled to non-zero angular momenta [31] are mainly responsible for the quenching of 
the Gamow-Teller strength in the ground states. However, we also observed, by performing 
calculations for two different residual interactions (the Brown-Richter interaction and the 
original Kuo-Brown interaction), that the calculated Gamow-Teller strength is rather sensi- 
tive to the residual interaction [13]. Moreover no systematic trends between the calculations 
using these forces and the data could be detected. As we demonstrate in the following, the 
situation improves significantly if one employs the modified Kuo-Brown interaction KB3. 

Our results for the total Gamow-Teller strengths B{GT^) = (((Tr_|_)^) are listed in Table 
1. (As our calculation obeys the Ikeda sum rule, values for 5(GT_) are readily obtained 
by adding 3{N — Z) to the B{GT^) values.) We observe that the calculated values are sys- 
tematically larger than the B[GT_^_) values deduced from intermediate-energy {n,p) charge- 
exchange cross section data at forward angles, which are known to be dominated by the GT+ 
operator and currently provide the only experimental determination of the Gamow-Teller 
strength function. However, such a systematic overestimation of the Gamow-Teller strength 
by shell model calculations is familiar from work in the srf-shell and can be attributed to 
an in-medium renormalization of the axial- vector coupling constant g^^ [3]. Since the (n,p) 
data are usually normalized to low-energy beta-decay rates, they are therefore indirectly 
also subject to any renormalization of qa- To account for this renormalization the shell 
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model results are usually multiplied by {q'a I QaY ~ = 0.64 [5,6]. If we apply this 

rescaling to our B{GT^) results for the /p-shell we find good agreement between the SMMC 
calculations and the data, as is shown in Fig. 4. For ^^Ti and "^^Ni our renormalized B{GT^) 
values deviate slightly from the measured Gamow- Teller strengths, indicating possible limits 
of the present model space. The slight discrepancy in the case of ^^Ti might be of some 
importance for shell model calculations of the double-beta decay rate of ^*^Ca [32]. 

^From Fig. 4 we conclude that i) Full Ohoj shell model calculations describe the system- 
atics of the Gamow- Teller quenching in the /p-shell; ii) Reproduction of the data requires a 
renormalization by 0.64, in agreement with the usual assumption of an in-medium modifica- 
tion of Qj^; and iii) these results for the /p-shell are consistent with those deduced previously 
for the sd-shell and indicate that both conclusions i) and ii) might be universal. If one 
accepts these conclusions, the agreement between data and theory suggests both that the 
(n,p) experiments do not miss any significant strength at higher energies and that the KB3 
interaction well-describes isoscalar correlations in /p-shell nuclei. 

The unrenormalized B[GT_^) strengths for the various isotopic chains are plotted in 
Fig. 5. As predicted by the simple single-particle estimate [33,34], the B{GT^) strength is 
roughly constant in the titanium isotopes, corresponding to a quenching factor of 3. For the 
other isotopic chains our calculation is in agreement with a recently suggested systematics 
of the experimental data [34] in which the total B{GT^) strength for mid- /p-shell nuclei is 
proportional to the numbers of valence protons and neutron holes in the /p-shell. As can 
be seen in Fig. 5, the linear proportionality to the number of neutron holes is also found in 
our Monte Carlo results. We also observe that the Gamow- Teller strength in the nuclei with 
12 neutron holes, corresponding to the magic neutron number = 28, is larger than the 
trend in the nuclei of the same isotope chain, which might further indicate an overestimate 
of the shell closure by the KB3 interaction. /,From the slopes of the B[GTj^) strengths for 
the various isotope chains we find that B{GT^) scales with the number of valence protons, 
as expected. 

For the even-even N = Z nucleus ^^^Ge we find an unrenormalized total Gamow-Teller 
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strength of B{GT+) = 7.91 ± 0.54. For the odd-odd N = Z nuclei we calculate B{GT+) 
values of 8.1 ±2.5 {^^Mn), 9.1 ±1.7 (^^Co), 6.6 ±2.8 (^^Cu) and 9.1 ±2.2 (^^Ga). The large 
uncertainties prohibit us from drawing any meaningful conclusions from these values. 



It is well known that pairing between like nucleons plays an essential role for the ground 
state properties of even-even nuclei. In a first approximation the pairing can be described 
by the BCS model, which assumes that like nucleons are coupled to J = pairs. We have 
studied the BCS-like pairing content of the ground states by measuring the expectation 



The results, obtained after subtraction of the independent particle model values (10), are 
shown in Fig. 6 for three different isotope chains. As expected, the excess pair correlations 
are generally quite strong and exceed the independent particle model values by factors of 2 
to 4, reflecting the known coherence in the ground states of even-even nuclei. It is interesting 
to note that the proton pairing fields are not constant within an isotope chain, but usually 
increase with neutron number, showing that there are important neutron-proton correlations 
present in these ground states. The shell closure at = 28 is manifest in the neutron pairing. 
For ^^'^"'Fe and ^^Ni the neutrons prefer to take advantage of closing the /7/2-subshell, making 
the excess neutron pairing rather small (about a factor of two compared to the independent 
particle model). But once there are extra neutrons outside the closed subshell, the excess 
of pairing increases drastically. As is demonstrated in Fig. 7 these strong changes are not 
present in the average occupation numbers rij. For example, the proton occupation numbers 
show little variation within the iron chain ^^^^^Fe, although the occupation of the /7/2 orbit 



E. Proton-proton and neutron-neutron pairing 



values for the pairing fields, (A^A), for proton-proton and neutron-neutron pairing. Here 
the pair operator is defined as A^ = X)j,m>o '^jm'^jmi where m is the time-reverse orbit of m. 
For a Fermi gas with occupations rij = YjmWjm'^im) we note that 




2 
n ; 



(10) 
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is largest in the semi-magic nucleus Fe. A similarly smooth behaviour is found in the 
neutron occupation numbers, where the additional neutrons are, on average, added to the 
lowest possible orbital. It is interesting to note that, in ^^'^^Fe, as a result of a strong 
/7/2-/5/2 coupling, the occupation of the energetically unfavored f^j2 orbital is larger than 
that of the orbitals. The presence of neutron pairs in the orbital, as in 56,58pg^ 
increases the occupation number for protons in this orbital, at the expense of the occupation 
of the ground state orbital. The comparison of Figs. 6 and 7 is yet another example for the 
importance of nucleon correlations in the ground state beyond the mean field level. 

IV. CONCLUSIONS 

Using the recently developed shell model Monte Carlo approach we have studied the gross 
properties of even-even nuclei in the mid-/p shell within full {)huj shell model calculations. 
Our studies use the KB3 interaction, a minimally corrected version of the original Kuo-Brown 
G-matrix interaction. Conventional diagonalization approaches have already established 
that shell model calculations using this KB3 interaction give a very satisfying description 
of nuclei at the beginning of the /p-shell [5,6]. The present calculations supplement these 
studies and show that shell model calculations, using this physical interaction, satisfyingly 
describe the gross properties of (even-even) nuclei throughout that part of the /p-shell [A < 
64) where the influence of the gg/2 orbital should still be small. We view this reproduction of 
such a large body of data, as presented in this paper, as a remarkable success of a microscopic 
model. 

We recall that performing full {Ohoj) calculations of mid- /p-shell nuclei, as discussed 
in this paper, was unthinkable only two years ago. Thus, the present studies and their 
results are also a successful confirmation of the Monte Carlo approach to the nuclear shell 
model, and establish it as a powerful tool with which to study gross nuclear properties. It is 
certainly a useful complement to the conventional diagonalization approaches, which remain 
the method of choice for detailed spectroscopic properties. 
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Our calculations have shown that shell model calculations with the KB3 interaction 
describe the binding energies of even-even and N = Z nuclei with = 48 — 64 to within 1 
MeV or better. We are quite aware that shell model calculations often do even better. It is 
conceivable that an improved reproduction of experimental masses requires the correction 
of tiny residual monopole defects in the KB3 interaction, as suggested in Ref. [5]. Such 
an optimization of the interaction is beyond the scope of the present calculation, but is 
conceivable with present computer capabilities. Our studies also indicate that the KB3 
interaction slightly overemphasizes the = 28 shell closure. 

Using conventional effective charges that account for coupling to the giant quadrupole 
resonance outside our model space, our calculated B[E2) strengths apparently reproduce 
the trend in the data suggested by the experimental B{E2) values for the transition from 
the ground state to the first 2+ state. This transition exhaust typically 70 — 80% of the total 
B(E2) strength in even-even mid-/p-shell nuclei. 

The most important result of our present study is a clarification of the Gamow-Teller 
quenching puzzle in the astrophysically important mid-/p-shell region, which has been the 
focus of attention for several years. While several mechanisms for the Gamow-Teller quench- 
ing have already been identified [30,31,13], we are, for the first time, able to demonstrate 
that the experimentally observed quenching is consistently reproduced by the correlations 
within the full /p-shell, if the standard renormalization factor of 0.64 is invoked. We find 
that the complete Ohuj calculations recover significantly more quenching than truncated 
shell model studies. Our results are consistent with previous (Ohuj) calculations of .sd-shell 
nuclei [3] which established that the Gamow-Teller strength is quenched by an additional 
renormalization factor 0.64 beyond the many-body correlations within that full model space. 
Thus, this additional renormalization of the Gamow-Teller strength appears to be universal 
and likely originates outside of nuclear configuration mixing within one major shell; it is 
consistent with = 1 in the nuclear medium. 

This work was supported by NSF grants PHY91-15574, PHY94-12818 and PHY94-20470 
and by DOE grant DE-FG02-91-ER40608. We thank Prof. A. Poves for useful discussions 

14 



and for supplying us with some results of his shell model diagonalizations prior to publication. 
Computational cycles were provided by the Concurrent Supercomputing Consortium and by 
the VPP500, a vector parallel processor at the RIKEN supercomputing facility; we thank 
Drs. I. Tanihata and S. Ohta for their assistance with the latter. 



15 



REFERENCES 

[1] O. Haxel, J.H.D. Jensen and H.E. Suess, Phys. Rev. 75, 1766 (1949). 
[2] S. Cohen and D. Kurath, Nucl. Phys. 73, 1 (1965). 

[3] B.A. Brown and B.H. Wildenthal, Ann. Rev. Nucl. Part. Sci. 38, 29 (1988) 
[4] T.T.S. Kuo and G.E. Brown, Nucl. Phys. A114, 241 (1968) 

[5] E. Caurier, A. P. Zuker, A. Poves and G. Martinez-Pinedo, Phys. Rev. C50, 225 (1994) 
[6] E. Caurier, A. Poves and A. P. Zuker, Phys. Rev. C, in print 

[7] N. Auerbach, G. F. Bertsch, B. A. Brown and L. Zhao, Nucl. Phys. A556, 190 (1993) 

[8] M. B. Aufderheide, S. D. Bloom, D. A. Rcslcr and G. J. Mathews, Phys. Rev. C47, 
2961 (1993) 

[9] A. L. Williams et al., submitted to Phys. Rev. C 
[10] A. Poves and A. P. Zuker, Phys. Rep. 70, 235 (1981) 

[11] G.H. Lang, C.W. Johnson, S.E. Koonin, and W.E. Ormand, Phys. Rev. C48, 1518 
(1993) 

[12] Y. Alhassid, D. J. Dean, S. E. Koonin, G. Lang and W. E. Ormand, Phys. Rev. Lett. 
72, 613 (1994) 

[13] D. J. Dean, P. B. Radha, K. Langanke, Y. Alhassid, S. E. Koonin and W. E. Ormand, 
Phys. Rev. Lett. 72, 4066 (1994) 

[14] D. J. Dean. S. E. Koonin, K. Langanke, B. P. Radha and Y. Alhassid, Phys. Rev. Lett., 
in print 

[15] C.W. Johnson, S.E. Koonin, G.H. Lang and W.E. Ormand, Phys. Rev. Lett. 69, 3157 
(1992) 

16 



[16] W.E. Ormand, D.J. Dean, C.W. Johnson, G.H. Lang and S.E. Koonin, Phys. Rev. C49, 

1422 (1994) 

[17] J. Hubbard, Phys. Rev. Lett. 3, 77 (1959); R.L. Stratonovich, Dokh Akad. Nauk. SSSR 
115, 1097 (1957) 

[18] A. Poves, private communication 

[19] Y. Alhassid, D.J. Dean, W.E. Ormand, S.E. Koonin and K. Langanke, unpublished. 
[20] A. Richter, Nucl. Phys. A507, 99c (1990) 

[21] D.L Sober, B.C. Metsch. W. Kniipfer, G. Eulenberg, G. Kiichler, A. Richter, E. Spamer 
and W. StefFen, Phys. Rev. C31, 2054 (1985) 

[22] W. Mettner, A. Ricliter, W. Stock, B.C. Metsch and A.G.M. Van Hees, NucL Phys. 
A473, 160 (1987) 

[23] A. Arima, K. Shimizu, W. Bentz and H. Hyuga, Advances in Nuclear Physics 18, eds. 
J.W. Negele and E. Vogt, (Plenum Press, New York, 1987) p. 1 

[24] W.P. Alford et al, Nucl. Phys. A514, 49 1990 

[25] M.C. Vetterli et ai, Phys. Rev. C40 40, 559 (1989) 

[26] S. El-Kateb et al, Phys. Rev. C49, 3129 (1994) 

[27] S. Raman et al., Atomic Data and Nuclear Data Tables 36, 1987 1 

[28] Nuclear Data Sheets (Academic Press, San Diego) 

[29] G. Kraus et al, Phys. Rev. Lett. 73, 1773 (1994) 

[30] J. Engel, P. Vogel and M. R. Zirnbauer, Phys. Rev. C37, 731 (1988) 

[31] M. Hino, K. Muto and T. Oda, Phys. rev. C37, 1328 (1988) 

[32] E. Caurier, A. Poves and A. P. Zuker, Phys. Lett. B252, 13 (1990) 

17 



[33] G. M. Fuller, W. A. Fowler and M. J. Newman, Ap. J. 252, 715 (1982) 
[34] S. E. Koonin and K. Langanke, Phys. Lett. B326, 5 (1994) 
[35] T. Ronnqvist et ai, Nucl. Phys. A563, 225 (1993) 



18 



TABLES 

TABLE L Total B{M1) (in B{E2) (in e^fm'') and unrenormalized B(GT+) strengths as 

calculated in the SMMC approach. For comparison the experimental B{GT^) strengths [24-26,9] 
and the B{E2) values for the 0'^ 2^ transition [27,29] are also given. 



nucleus 


5(M1) 


B{E2){SMMC) 


B(E2) (exp) 


B{GT+) (SMMC) 


BiGT+) (exp) 


48Ti 


10.2 


± 


1.2 


455 ± 25 


720 ± 40 


1.13 


± 


0.18 


1.31±0.2 


50Ti 


12.5 




1.0 


415 ± 50 


290 ± 40 


1.47 


± 


0.16 






12.5 


± 


1.0 


465 ± 55 


> 250 


1.11 


± 


0.16 






13.5 


± 


1.5 


450 ± 80 




0.97 


± 


0.21 




48Cr 


13.8 


± 


1.7 


945 ± 45 


1330±200 


4.37 


± 


0.35 




50Cr 


14.5 


± 


2.5 


890 ±90 


1080±60 


3.51 


± 


0.27 




52Cr 


18.9 




2.2 


645 ±75 


660±30 


3.51 


± 


0.19 






13.0 


± 


2.5 


890 ±90 


870±40 


2.21 


± 


0.22 






16.2 


± 


2.0 


840 ±90 




1.50 


± 


0.21 




52Fe 


18.9 


± 


1.4 


1055 ± 50 




7.10 


± 


0.42 






16.5 


± 


2.8 


750 ± 80 


620±50 


6.05 


± 


0.45 


3.1±0.6 




20.4 




3.0 


990 ± 65 


980±40 


3.99 


± 


0.27 


2.85±0.3 


S^Fe 


20.3 


± 


3.0 


1010 ± 65 


1200±40 


3.06 


± 


0.28 






17.3 


± 


2.1 


1105 ± 65 


930±180 


1.80 


± 


0.24 




56^1 


22.5 


± 


1.2 


515 ± 40 


600±120 


9.86 


± 


0.38 




58^1 


20.0 


± 


2.0 


960 ± 75 


695±20 


6.72 


± 


0.50 


3.76±0.4 


60Ni 


22.0 




2.5 


1065 ± 75 


935±15 


5.18 


± 


0.39 


3.11±0.08 


62Ni 


19.6 


± 


2.9 


1010 ± 85 


890±25 


3.43 


± 


0.40 


2.53±0.07 




18.9 


± 


2.7 


1165 ± 75 


760±80 


1.73 


± 


0.29 


1.72±0.09 


60Zn 


19.5 


± 


1.2 


1335 ± 50 




8.13 


± 


0.39 




62Zn 


19.0 




2.2 


1350 ± 70 


1230±90 


5.60 


± 


0.50 




64Zn 


23.6 




2.2 


1225 ± 65 


1440±120 


4.13 


± 


0.34 
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FIGURES 

FIG. 1. Comparison of the total Gamow- Teller strengths B{GT^) for ^^'^^Fe in a series of 
direct diagonalizations with decreasing level of truncation [6] with present full /|j-shell result (solid 
symbols plotted at t = 10). The open symbols at t = 10 show the extrapolated no-truncation 
result of [6]. 

FIG. 2. Upper panel (a): Comparison of the mass excesses AM as calculated within the SMMC 
approach with the data. Lower panel (b): Discrepancy between the SMMC results for the mass 
excesses and the data, 5AM. The solid line shows the average discrepancy, 450 keV, while the 
dashed lines show the rms variation about this value. 

FIG. 3. Comparison of experimental B(E2, 0^ — >• 2^) strengths with the total B(E2) strength 
calculated in the SMMC approach for various fp-shell nuclei with either proton or neutron number 
N = 28. 

FIG. 4. Comparison of the renormalized total Gamow- Teller strength, as calculated within 
the present SMMC approach, and the experimental B(GT^) values, deduced from {n,p) data 
[24-26,9,35]. Note that the two measurements of B(GT+) for ^''Fe summed the strength up to 8 
MeV (3.1 ± 0.6, [25]) and up to 9 MeV (3.5 ± 0.7, [35]). 

FIG. 5. Unrenormalized total Gamow- Teller strength for various isotope chains as a function 
of the number of neutron holes (TV — 20) in the fp-shell. 

FIG. 6. Expectation values for the proton (upper) and neutron (lower) pairing fields, as calcu- 
lated in the even iron, nickel and zinc isotopes. The values of (A^A) for the independent particle 
model have been subtracted. 

FIG. 7. Occupation numbers of the various orbitals in the iron isotopes ^^^^^Fe, as calculated 
in the present SMMC approach. For clarity the proton occupation numbers of the p^f2i /5/21 
Pl/2 orbitals have been multiplied by a factor of 5. 
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